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ABSTRACT 

I provide estimates of the ultraviolet and visible light luminosity density at 
z ~ 6 after accounting for the contribution from faint galaxies below the detection 
limit of deep Hubble and Spitzer surveys. I find the rest-frame V— band luminosity 
density is a factor of ~ 2 — 3 below the ultraviolet luminosity density at z ~6. 
This implies that the maximal age of the stellar population at z ~ 6, for a 
Salpeter initial mass function, and a single, passively evolving burst, must be 
<100 Myr. If the stars in z ~ 6 galaxies are remnants of the star-formation 
that was responsible for ionizing the intergalactic medium, reionization must 
have been a brief process that was completed at z < 7. This assumes the most 
current estimates of the clumping factor and escape fraction and a Salpeter slope 
extending up to 200 M for the stellar initial mass function (IMF; dN/dMocM a , 
a = —2.3). Unless the ratio of the clumping factor to escape fraction is less than 
60, a Salpeter slope for the stellar IMF and reionization redshift higher than 7 
is ruled out. In order to maintain an ionized intergalactic medium from redshift 
9 onwards, the stellar IMF must have a slope of a = —1.65 even if stars as 
massive as ~200 M are formed. Correspondingly, if the intergalactic medium 
was ionized from redshift 11 onwards, the IMF must have a ~ —1.5. The range 
of stellar mass densities at z ~ 6 straddled by IMFs which result in reionization 
at z > 7 is 1.3±0.4xl0 7 M Mpc" 3 . 

Subject headings: galaxies: stellar content - - galaxies: high-redshift — early 
universe 



1. Introduction 

Deep optical/infrared imaging surveys with Hubble and Spitzer are providing an un- 
precedented opportunity to study star-forming galaxies at high redshifts. It is increasingly 
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clear that these star-forming galaxies are responsibl e for reionizing the intergalactic m edium 
(IGM), a process which seems to complete by z ~ 6 (IBecker et al.ll200ll ; iFan et al.ll2006l ). The 
shape of the luminosity function of high redshift active galactic nuclei detected by the Sloan 
Digital Sky Survey (SDSS) indicates that ionizing photons from supermassi ye black holes fa ll 
many orders of magnitude short of the minimum required to ionize the IGM (IFan et al.ll200ll ). 
Even the curre n t star - forma tion in Lyman-break galaxies (LBGs) at z ~ 6 as measured by 
Bouwens et al.l (120061 . 120071 ). falls about a factor of 6 — 9 below the minimum required to 
maintain an ionized IGM for canonical estimates of two key parameters, the clumping factor 
of the gas and the escape fraction, assuming a Salpeter initial mass function (IMF). One 
possibility is that the contribution from star-formation in faint galaxies, below the detection 
threshold of current surveys, is higher than has been estimated. This translates to a faint 
end slope for the galaxy ultraviolet luminosity function of —1.9, a value that is conside rably 
steeper than the value of —1.74 that has currently been measured (IBouwens et al.ll2007l ). Ev- 
idence for such a steep slope is present in estimates of star-formation rate densities derived 
from the comoving gamma-ray burst number density (jChary et al.ll2007al ). However, that 
measurement has large uncertainties limited by the small number of high redshift GRBs cur- 
rently known and the uncertain conversion between GRB number density and star-formation 
rate density. 

An alternate possibility is that the star-formation rate in galaxies currently detected at 
z ~ 6 was higher in the past. This implies that the ultraviolet luminosity function of galaxies 
must show luminosity evolution between z ~ 6 and higher redshifts. Although constraints 
on the ultraviolet luminosity function at higher redshift are difficult to obtain from existing 
data, the upp er limits on bright L yman-break galaxies at z » 6, seem to rule out such 
an evolution (IBouwens et al.l 120041 ) . In fact, the upper limits seem to suggest a negative 
evolution with increasing redshift, in the sense that galaxies at the bright end of the UV 
luminosity function at z ~ 6 were fainter at higher redshift in the ultraviolet. There are 
no available constraints on the evolution of the faint end of the UV luminosity function at 
z > 6. 



Measurement of the stellar mass density provides an alternate constraint on the past 
history of star- formation. Deep Spitzer surveys which detect the rest-frame visible light from 
high r edshift galaxies a l low the stellar ma ss in galaxies at z ~ 6 to be measured (jYan et al. 



20061 ; IStark et al.l 12006c lEyles et al.l 120051 ) . The galaxies which are detected in rest-frame 
optical light, are primarily the brightest galaxies in the ultraviolet. The stellar mass density 
directly inferred from the Spitzer detections and stacking analysis are therefore, a lower limit. 
By making reasonable assumptions for the rest-frame V— band luminosity of galaxies which 
are undetected in the Spitzer data but seen in the ultraviolet, it is possible to estimate the 
true co-moving stellar mass density at z ~ 6. Since the stars at z ~ 6 are remnants of 
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star-formation at higher redshifts, we can use this stellar mass density estimate to constrain 
the number of ionizing photons produced by star-formation in the past for comparison with 
the minimum number of photons required to keep the IGM ionized. 

At present, the exact redshift of reionization is not known. Wilkinson Microwave 
Anisotropy Probe (WMAP) polarization measurements indicate a large optical depth to scat- 
tering of microwave background photons suggesting that the intergalactic medium (IGM) 



was i onized at some time between redshifts of 7 and 14 (IKogut et al.l 120031 ; ISpergel et al. 



2003). 



In this paper, we estimate the optical and ultraviolet luminosity density at z ~ 6 after 
correcting for faint galaxies undetected by current surveys. We then calculate the stellar 
mass density and average age of the stellar population in galaxies, required to reproduce 
the luminosity density estimates for different shapes of the stellar IMF. We also quantify 
the number of ionizing photons produced over the past history of the starburst and assess 
if they are sufficient to ionize the IGM at z > 6. Finally, we constrain the stellar initial 
mass function at z > 6 as a function of the reionization history of the Universe. A standard 
f2 M =0.27, ^a=0.73, H =71 km s -1 Mpc -1 cosmology is adopted throughout this paper. 



2. Optical and Ultraviolet Luminosity Density at z~6 



Deep Spitzer imaging data which detects the rest-frame visible emi ssion of z > 6 gal axies 
has yielded the fi rst constraints on the stellar mass density at z ~ 6. lYan et al.l (120061 ) and 
Eyles et al.l (120051 ) measured the IRAC photometry of Lyman-break galaxies selected using 
the i— dropout technique in the Great Observatories Origins Deep Survey (GOODS) fields. 
By fitting the multiband photometry of each galaxy with population synthesis models and 
deriving the stellar mass of each galaxy, they find a lower limit to the stellar mass density of 
0.1 — 0.7 x 10 7 M Mpc -3 , assuming a Salpeter initial mass function. The value is a lower 
limit since only the brightest galaxies in the ultraviolet are typically detected while many 
others ar e blended w i th bri ghter foreground sources in the IRAC images. At lower redshifts 
of z ~ 5, IStark et al.l (120061 ) estimate the stellar mass density to be 0.6 — 1 x 10 7 M^ Mpc ~ 3 



by com bining a spectroscopic and photometric redshift sample in GOODS-S. IChary et al. 



( 120071 ) adopt an alternate technique where they fit for the stellar masses of only a spectro- 
scopic sample of galaxies at 5 < z < 6.5 to derive a lower limit to the stellar mass density 
of 0.3— l.lxlO 6 Mfn Mpc~ 3 at z ~ 6. They adopt the ultraviolet luminosity function of 



Bouwens et al.l ( 120071 ) and make the assumption that the rest-frame UV/V— band luminosity 
ratios of the Spitzer undetected galaxies at these redshifts are the same as the detected galax- 
ies. The Spitzer detected galaxies have z/L^(1500A)> 2xl0 10 L and z/L^(5500A)> 10 10 L . 
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The median visible to ultraviolet luminosity ratio for the Spitzer detected galaxies is 0.5. 
They then apply a completeness correction by scaling the UV luminosity density of their 
spectroscopic sa mple to the total UV luminosity density of Lyman-break galaxies at z ~ 6 
as estimated by lBouwens et al.l (120071 ). This completeness correction is a factor of 34 which 



results in a total stellar mass density of 1 — 3.7 x 10 7 M Mpc -3 at z ~ 6 for a Salpeter IMF. 
The corresponding UV and V-band luminosity densities derived by integrating down to a UV 
luminosity of 10 7 L is 1.4xl0 8 L Mpc -3 and 6.6xl0 7 L Mpc~ 3 , respectively. We call this 
the "high-V" case. Integrating to zero luminosity results in only a 10% upward correction to 
these values. These values are in excellent a greement wi t h the stellar mass density at z ~ 6 
obtained in the hydrodynamic simulations of iDave et al.l (12006 ) but substantially lower than 



the value of 8.8xlO 7 M Mpc -3 derived bv iNagamine et al. I (120081 ) 



The assumption that the visible to ultraviolet luminosity ratios of faint LBGs is similar 
to that of brighter LBGs is subject to some uncertainty. Faint LBGs are predominantly 
in low mass halos. The supernova feedback from the initial burst of star-formation might 
be sufficient to inhibit or delay further star-formation. This would suggest that UV faint 
galaxies have lower visible to ultraviolet luminosity ratios than galaxies at the bright end of 
the UV luminosity function. This trend between rest-frame optical -ultraviolet colors with 



ultraviolet luminosity appears to be seen in semi-analytical models (jSomerville et al.l 12001 
personal communication). To assess the dependence of the z ~ 6 luminosity density on 
the optical luminosities of faint galaxies, we assume that the UV/V— band luminosity ratios 
of faint Spitzer undetected galaxies follow the correponding ratios seen in semi-analytical 
models. The ratio of the visible to ultraviolet luminosity density in the models, for galaxies 
below the Spitzer detection threshold, is 0.28. For the UV bright galaxies, some of which 
are undetected due to blending, we adopt the UV/V— band luminosity ratios of the Spitzer 
detected sources which are similar to the colors of UV bright galaxies in the semi-analytical 
models. We find that the resultant V-band luminosity density under these assumptions, 
obtained by integrating dow n to ;/L, y (1350A)~ 10 7 L is 4.2 xlO 7 L Mpc -3 which is 30% 
smaller than the estimate of IChary et al.l (120071 ). We call this the "low-V" case. 



More recently, massive, evolved stellar populations have been claimed in candidate high 
redshift galaxies which would boost the derived V — band luminosity densities and therefor e 
the stellar mass densities, by factors of several (IWiklind et al. 1120081 ; iPanagia et al. 1 120051 ) . 
However, since a significant fraction of such sources are detected at 24/zm, a wavelength which 
traces the dust em ission from galaxies , the evidence is strongly against these objects being 
at high redshifts ( Chary et al.l l2007bl )n. Furthermore, through mid- infrared photometric 



lr The exception is a galaxy at z=5.554 which might have its stellar mass estimate erroneously elevated 
due to the presence of Ha in the Spitzer passbands. See source ID44 in Chary et al. (2007) for details. 
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redshifts, it has been demonstrated that the best candidate am ong these proposed high 



redshift galaxies is a dusty, infrared luminous galaxies at z ~ 1.7 (IChary et al.l l2007bl ) . We 
ignore this population of galaxies for the rest of the discussion until there is any reliable 
evidence that they really are at z » 5. 

Typically, stellar mass densities are estimated by fitting the multi-wavelength spectral 
energy distribution of individual galaxies with population synthesis models and integrating 
over the resultant stellar mass distribution. However, the stellar mass density is sensitive 
to the choice of the stellar mass function, which is one of the parameters we are trying to 
constrain. Furthermore, since we do not have multiwavelength detections of faint galaxies, 
the only available approach involves fitting a population synthesis model to the integrated 
rest-frame UV and V-band luminosity density estimates at z ~ 6 derived above. 

We adopt the Starburst99 set of models which have the Pado va tracks for the ther- 



mally pulsing asymptotic gi ant branch stars (TP-AGB) incorporated (ILeitherer et al. 111999 



Vazquez fc Leitherer II2005I ). The variables in the models are the slope of the stellar IMF, 
metallicity, age of the stellar population and stellar mass density. We consider both O.O2Z 
and 0AZ Q templates and single burst models. We also con volve the templates with a model 



for the absorption by the Ly— a forest (IMadau et al. Ill995l ). The fit to the luminosity density 



constrains the stellar mass density as well as the total number of ionizing photons produced 
over the duration of the starburst. Adopting a single burst model yields maximal ages and 
stellar mass densities while the total number of ionizing photons produced is insensitive to 
the time history of star-formation. 

If we adopt a Salpeter IMF with metallicity 0.02 Z , we find that the optical and ultra- 
violet luminosity density at z ~ 6 are fit by a stellar population of age 100 Myr and with 
a stellar mass density of 1.6xl0 7 M Q Mpc~ 3 (Figure 1). If instead we consider the "low-V" 
case, we derive a maximal stellar age of 50 Myr and a total stellar mass of 0.8 x 10 7 M Mpc~ 3 . 
The presence of extinction at z ~ 6 would imply that reproducing the observed optical to 
ultraviolet luminosity densities will require a bluer intrinsic stellar population, corresponding 
to a younger age. Since the average extinction in star-forming regions at z ~ 6 is unknown, 
we do not attempt to obtain fits that include extinction. 

Changing the stellar initial mass function will change the derived stellar mass densities. 
Since we use the optical and ultraviolet luminosity densities at z ~ 6 as the overall constraint, 
and evolve a single star burst of fixed mass, the true stellar mass density is the normalization 
factor for the fixed mass burst. Thus, for every stellar IMF that we consider, we can derive 
the age, stellar mass density and number of ionizing photons produced by requiring that the 
luminosity density estimates at z ~ 6 are reproduced. We note that the age is not a strongly 
constrained parameter since we are fitting to only observations at two wavelengths and the 
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exact age would strongly depend on the average time history of star-formation. By adopting 
a single burst model, we only derive maximal ages for a single stellar population. 

The ultraviolet and optical luminosity density at z ~ 6 are dominated by ongoing star- 
formation and intermediate mass stars respectively, the latter depending on the age of the 
starburst. For a particular IMF, the stars with lifetimes shorter than the age of the starburst 
at that redshift would have presumably evolved to their end states i.e. neutron stars or black 
holes. They would still be contributing to the total stellar mass density in the models as 
well as the total number of ionizing photons produced by the starburst, even if they dont 
contribute to the ultraviolet/optical luminosities at z ~ 6. 

An upper limit to the stellar mass density for each IMF can also be derived by adopting a 
two component stellar population synthesis model (Figure 1, right panel). The young stellar 
population, which dominates the UV luminosity density is assumed to be ~10 Myr old. An 
old, high mass-to-light ratio stellar population is assumed to form in a single burst at z ~ 20 
and evolve passively to z ~ 6. This implies a stellar age of 770 Myr. Naturally, a stellar 
IMF in which all stars have lifetimes younger than 770 Myr will yield unphysical stellar 
masses and is represented with an ellipsis i n Table 2. We not e that such a two component 



starburst has been considered previously by I Chary et al.l (120071 ) who showed that the number 



of ionizing photons produced even in such an extreme scenario falls short of maintaining an 
ionized IGM, assuming a Salpeter IMF. Although concentrating all the star-formation into 
a single z ~ 20 burst increases the total number of ionizing photons produced by factors of 
a few (see Table 1 and 2), the recombination rate increases as (1 + z) 3 which implies that a 
much larger number of ionizing photons are required to keep the Universe ionized at these 
redshifts. 



3. Constraining Reionization 

By fitting the measured optical and ultraviolet luminosity density at z ~ 6 with a 
evolving, single burst model, we can determine the age of the stellar population and the 
stellar mass density. Integration of the number of photons produced shortward of the Lyman 
limit over the age of the population, yields the total number of ionizing photons produced 
by a starburst template which fits the luminosity density values. The integral constraint 
has the advantage of being independent of the exact star-formation history which is likely 
to have a multitude of ages and e-folding times. Since the vast majority of ionizing photons 
are produced in massive stars, the single burst model also has the advantage that the time 
integral of the number of ionizing photons, for starburst ages larger than ~10 Myr, converges 
to a single value. 
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The number of ionizing photons produced over the history of star-formation can be 
compared with the number of ionizing photons required to ionize the IGM. The number of 
photons required to initiate reionization is only 1 photon per baryon. However, to account 
for recombination, the production rate of ionizing photons per baryon must be ~10 times 
higher. To quantify the rate of ionizing photons N(z) required to keep the IGM ionized at 
any redshift, the rate of ionizing photons must be greater than or equal to the recombination 
rate. 

The recombination rate R is: 

R = n e nun ®b C s -1 Mpc -3 (1) 

where n e is the comoving electron density, tihii is the comoving ionized hydrogen density, 
is the recombination coefficient to excited states of hydrogen while C is the ionized hydrogen 
clumping factor defined a s < n^n >/< ^htt > 2 - Since Hell reionization takes place at much 



lower redshifts z ~ 3 — 4 (ISokasian et al 



2002 ). the comoving electron density is: 



n e = n m i + n He ii (2) 

and not n H n+2nH e as has previously been adopted, although the difference is negligibly 
small because of the small fractional number density of He. We adopt a j B=2.59xlO~ 13 (l + 
z) 3 cm 3 s _1 for the comoving recombination coefficient (Case B), assuming a gas temperature 
of 10 4 K. So, the photon production rate per baryon to balance the recombination rate at 
any redshift, is: 

N(z) = - s- 1 (3) 
n h 

The number of photons per baryon to start reionization at any redshift N (z) is simply the 
ionized hydrogen fraction at that redshift. 

JVoW = = (4) 

So, the total number of ionizing photons per baryon required to maintain the increasing 
ionized fraction of the IGM between any redshift z and redshift z, where Zq > z, is: 

r z ' =z ■ dt 

N(z) = N (z) + J _ N(z')—dz' (5) 

Reionization of hydrogen is complete by z = 6 while reionization of helium is not. So, 
Nq(z = 6) is the fractional volume density of hydrogen which is 0.93. Equation 5 can be 
rewritten as: 

r z ' =t > dt 
N{z = 6) = 0.93 + 6.475 x lO" 20 / C X mi (xhii + XhcIi) (1 + z'f —dz' (6) 

Jz'=z dz 
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Since VL^h 2 = 0.0223, the comoving baryon volume density is 0.25 baryons m -3 . xnn and XHeii 
are the volume averaged gas fractions for ionized hydrogen and ionized helium respectively. 
Care should be taken that these are derived with respect to the number density of baryons. 
Typical simulations provide ratios with respect to the number density of neutral hydrogen 
or neutral helium. In that case, we have adopted nn=0.9268nb and nH e =0.073nb. 

If we neglect the small contribution to the electron density from Helium, Equation 6 
can be simplified to: 



N(z = 6) « 0.93 + 6.475 x 10" 



20 



I X—ZQ 



'C {l + z'f ^-dz' 
dz' 



(7) 



where C is < Wjj n 



> / < rib > 2 , a quantity directly available from the llrac fc Cen I (120071 ) 
simulations (See their Figure 12). However, their estimate of C is an average over the entire 
comoving volume. In practise, the dumpiness of gas which is within the source-forming halos 
is already factored into the estimate of escape fraction. So, the clumping factor should be a 
measure of the dumpiness of the IGM excluding halos with star-formation in them (Figure 
2a). Therefore, we use the exact relation as shown in Equation 6. 

We adopt the most current estimate for the reionization history, redshift dependent 
clumping fac t or (F igure 2a), n e , nmi and nneii from the high resolution simulations of 
Trac fc Cen I (120071 ). A comparison between the number of photons per bar yon from these 



equat ions compared with a similar relation derived for complete reionization by lMadau et al 



(119991 ) is shown in Figure 2b. The number of ionizing photons required is sensitive to the 
reionization history since the recombination rate is proportional to the square of the density. 
If reionization is completed earlier, the recombination rate is proportionately higher than 
for a partially ionized medium. In that case, the number of photons per baryon needs to be 
higher to maintain the ionized state of the IGM as shown by the triple dot dashed line in 
Figure 2b. 



We note that the reionization history of iTrac &: Cen I (120071 ) corresponds to the "late" 
reionization history which would result in a W MAP optical depth t o Thomson scattering 
close to the measured value of r ~ 0.09 ± 0.03 JSpergel et al.~ 2007 ). As characterized in 



Greif fc Bromm ( 2006 ) and Wyithe fc Loeb ( 20031 ) . shifting the process of reionization to 
z ~ 10 would result in optical depth values which are at the high end of the observed range. 

For each stellar initial mass function that we adopt, we obtain a library of template 
SED which are a function of the age of the starburst. We first determine the combination 
of age and mass which results in a best fit to the optical and ultraviolet luminosity density. 
This yields the stellar mass density at z = 6 for that IMF. The range of age and mass values 
for a subset of IMFs considered here is shown in Table 1. As emphasized earlier this is a 
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maximal age for a single stellar population since we adopt a single instantaneous burst. 



The Starburst99 model yields the number of hydrogen ionizing photons produced as a 
function of age of the starburst. The escape fraction is a measure of the number of ionizing 
photons escaping from star-forming r egions in galaxies, i nto the IGM. The esc ape fraction has 
been constrained to be ~ 10 — 15% (ISiana et al. 1 120071 ; IShapley et al. 1120061 ) in the Lyman- 
break galaxy population. Thus, we adopt escape fractions of 0.1 and 0.15. By integrating the 
number of ionizing photons over the derived stellar age, we can estimate the total number 
of photons produced over the lifetime of the evolving galaxies, multiply that with the escape 
fraction and compare the result with the number of ionizing photons required to keep the 
IGM ionized as determined from Equation 6. 



3.1. Salpeter IMF 

Figure 3a shows the total number of ionizing photons produced by a starburst for a 
Salpeter IMF which is a power-law of the form dN/dMocM Q , a = —2.3. The metallicity was 
0.02 Z , escape fraction of 0.1 and clumping factor shown by the solid black line in Figure 
2a. The calculations corresponds to the "high-V" case Different lower mass cutoffs were 
adopted for the Salpeter IMF as shown by the different colored lines while the upper mass 
cutoff was fixed at 200 M Q . The black line is the number of ionizing photons required to 
keep the IGM ionized between z < z' < 6. The point where the colored lines which denote 
the number of ionizing photons escaping from the star-forming regions cross the black line 
is the epoch beyond which reionization was complete. 

For a Salpeter IMF extending down to 0.1 M i.e. with no low mass cutoff, the IGM was 
ionized for only 60 Myr prior to z ~ 6. This is an unphysically small age and indicates that 
if the IMF has a Salpeter shape in z ~ 6 galaxies, reionization must have been an extremely 
inhomogeneous process. 

As the low mass cutoff increases from 0.1 to 10 M Q , the optical luminosity density 
produced for a fixed age for the starburst decreases. So, as the low mass cutoff increases, 
the age of the starburst needs to become younger to reproduce both the UV and optical 
luminosity density. As a result, the number of ionizing photons produced over the duration 
of a starburst by a Salpeter IMF with a mass cutoff at 10 M is lower than that with a mass 
cutoff at O.1M . 

Even in the case where the low mass cutoff is at 5 M , which produces the largest 
number of ionizing photons, the process of reionization has to be shorter than 140 Myr. 
This indicates that unless the ratio of clumping factor to escape fraction is significantly 
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lower than what the simulations seem to indicate, a Salpeter IMF and z re i on > 7 is ruled out. 

The "low-V" case, since it corresponds to a bluer, stellar population which has a younger 
age, would result in fewer ionizing photons over the lifetime of the starburst, making the 
conclusions derived above even more stringent. Thus, unless explicitly stated otherwise, our 
calculations below are done for the "high-V" case. 

3.2. Modifications to the High Redshift Stellar IMF 

Motivated by the WMAP results which suggest the possibility of reionization as early 
as z reion ~ 11, we investigate a range of initial mass functions which would produce sufficient 
ionizing photons to keep the IGM ionized at z » 7. 

As shown in Figure 3a, changing the lower mass cutoff of a Salpeter IMF cannot increase 
Zreion above redshift of 7. It appears necessary to change the value of the IMF slope a. 
The flatter the IMF slope, the larger the contribution of massive stars. Since most of the 
ionizing photons are produced by massive stars, even a small increase in the slope results 
in a significant increase in the number of ionizing photons produced. Changing the low 
mass cutoff does not change the number of ionizing photons substantially except for the 
renormalization by the total mass in the IMF. If a = —2, the IMF has equal mass in 
logarithmic mass bins. So for a value of a > — 2, z re i on is relatively insensitive to a change 
in the low mass cutoff in the IMF. 

Figure 3b shows the number of ionizing photons produced by an IMF which extends 
from 1 — 200 M and has a range of slopes a. For a = —2.3, z re i on falls between 6 < z < 7 
while for a = —2.0, sufficient ionizing photons are produced to initiate reionization at z ~ 11 
and complete the reionization process by z ~ 6. By comparing Figure 3b and Figure 4, it 
can be seen that z reion for a > — 2 is insensitive to the low mass cutoff of the IMF. The 
resultant stellar mass densities at z ~ 6 for the IMFs considered here are shown in Table 1. 

If the formation of high mass stars is suppressed, the number of ionizing photons pro- 
duced per baryon are reduced. If the cutoff is at 50 M instead of 200 M , the net effect 
is to increase the value of a. For example, a stellar IMF with a mass range of 1-50 M , 
a = —1.7 produces a similar number of ionizing photons as an IMF extending between 1-200 
M , a = -2.1. 

The choice of metallicity in the models is also an important parameter. For the cal- 
culations we present, the results are based on the 0.02 Z metallicity models which are the 
lowest available metallicity in Starburst99, although we also evaluated the results for the 
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0.4 Z metallicity templates as illustrated in Figure 3b. Lower metallicity models produce 
~ 2—3 times as many ionizing photons as higher metallicity models (Figure 3b). However, we 
note that the number of ionizing photons produced for zer o metallicity single burst templates 
which have normal IMFs as presented by ISchaerer I (120031 ). are very similar to the number of 
ionizing photons produced for 0.02 Z metallicity templates to within ~10-20%. This implies 
that a Salpeter IMF even with lower metallicity templates than our adopted ones cannot 
produce sufficient ionizing photons. Thus, although the exact slope of the derived IMF is 
sensitive to the metallicity of the gas collapsing into the stars, the fact that the metallicity 
enrichment from massive stars is a relatively quick process and the metallicity of GRB hosts 
at z > 5 which are presumably star-forming galaxies, appear to be ~ 0.01 — 0.1 Z , leads us 
to conclude that the choice of 0.02 Z metallicity templates is a reasonable approximation. 



3.3. Dependence between Stellar IMF and z re i on 

We have generated template SEDs for a range of stellar initial mass function slopes 
and stellar mass ranges which reproduce the observed optical and ultraviolet luminosity 
density at z ~ 6. We consider two cases, one in which the stellar IMF extends between 
1 — 200 M and a second in which the mass range is 1 — 50 M . As emphasized earlier, the 
choice of cutoff at the low mass end is insensitive to the result once a > —2. We assume 
that the metallicity is 0.02 Z in either case. We consider a range of — 1 > a > —2.6 and 
calculate the number of ionizing photons produced for each scenario. We then calculate 
the redshift at which the number of ionizing photons from the stellar population exceed the 
minimum required to maintain the ionized state of hydrogen. It is important to note that 
this is not simply a comparison with the solid black lines in Figure 3 and 4 which assume 
the particula r reion ization history (and thereby clumping factor) from the simulations of 



Trac fc Cen I (120071 ). 



If reionization were complete at higher redshifts, the recombination rate between the 
end of reionization and z ~ 6 would be higher because of the larger values of n e and tihtt- To 



estima te this, we simply shift the reionization history and clumping factors of iTrac fc Cen 



(120071 ) to higher redshifts and assume that the clumping factors remained constant once the 
IGM was completely ionized. The number of ionizing photons per baryon in the scenario 
where reionization is complete by z ~ 11 is plotted as the triple dot dashed line in Figure 
2. We note t hat this is a f actor of ~4 higher than the "late" reionization history which 



the models of ITrac fc Cen I (120071 ) provide. Furthermore, even if the clumping factor were 
to increase after the IGM was completely reionized, this would only increase the minimum 
number of photons per baryon required. Thus, our estimates are, at worst, a lower limit to 
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the number of photons per baryon required to maintain the ionized state of the IGM. 

We find that as the redshift of reionization increases, the value of a needs to increase 
(Figure 5). This is not surprising since a higher z ieion implies a larger number of ionizing 
photons which are preferentially produced in more massive stars. We also find that as the 
high mass end of the stellar IMF is cut, from 200 M to 50 M , the slope a needs to increase 
such that a larger fraction of massive stars make up the shortfall in ionizing photons. 

The relation between z re i on and IMF slope is also sensitive to the ratio of the clumping 
factor to the escape fraction. Reduction in the clumping factor decreases the number of 
ionizing photons required since the recombination rate is reduced. Similarly, a large escape 
fraction implies that a larger number of photons escape star-forming regions. So, a fewer 
number of ionizing photons need to be produced to reionize the IGM. As a result, in Figure 
5, we plot the relation between z reion and a for different ratios of C/f esc . A redshift dependent 
clumping factor such as that shown in Figure 2a and an escape fraction of 0.1 is shown as 
the green lines. The red lines show the estimates for a clumping factor of 30 and an escape 
fraction of 0.1. The blue line show the estimates for a clumping factor of 10 and an escape 
fraction of 0.15 while the cyan line is an extreme case with a clumping factor of 10 and escape 
fraction of 0.5. A similar comparison between the redshift when reionization is complete and 
the IMF slope is shown in Figure 6 for the "low-V" case. 

As can be seen for the "high-V" case, with a Salpeter IMF (a = —2.3), unless the ratio 
of clumping factor to escape fraction is much smaller than 60, the stars in z ~ 6 galaxies 
cannot be responsible for reionization at z > 7. For our best estimate of the clumping factor 
and escape fraction, if z re i on = 9, the slope a = —1.65 while if z reion = 11, a = —1.5. In either 
scenario, only 0.2% of all baryons at z ~ 6 are processed through stars and are responsible 
for ionizing the IGM. 



Our conclusions are therefore broadly consistent with the results of iTumlinson et al 



(120041 ) who argue against the need for a very massive star population at high redshift. 
Their arguments are based on the metal abundance ratios in metal-poor Galactic halo stars. 
However, we disagree with them in the sense that a Salpeter IMF even with a low mass 
cutoff, is inadequate for producing sufficient ionizing photons to ionize the IGM while fitting 
the observed ultraviolet and optical luminosity density at z ~ 6 (Figure 3). 

Such a low mass tu rno ver in the stellar IMF may indeed be present at z < 2 as has been 



argued by lDave I (120081 ) and Ivan Dokkum I (120081 ). The former attempted to reconcile specific 



star-formation rate estimates between models and observations while the latter fit the color 
evolution of massive early- type cluster galaxies at z < 1. We att e mpte d to fit the observed 



luminosity densities at z ~ 6 with the IMF proposed by iDave I (120081 ). We assumed that 
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the (1 + z) 2 evolution in the mass turnover of the stellar IMF extends up to z ~ 4, similar 
to the evolution proposed by Ivan Dokkum I (120081 ). Thus, the IMF was assumed to have a 
slope of -1.3 at 0.1 <M< 12 M and -2.3 at 200>M> 12 M . For the "high-V" case, I find 
that this IMF produces 3.25 ionizing photons per baryon, an age of the stellar population 
of 90 Myr and a stellar mass of 1.2xl0 7 M Q Mpc~ 3 . Unlike the Salpeter IMFs with a sharp 
cutoff shown in Fig ure 3, such an IMF can account for the the "late" reionization history of 



Trac fc Cen I (120071 ) and would maintain an ionized IGM at z < 7. However, for the "low-V" 



case, this IMF would produce only 1.0 ionizing photons per baryon, result in a stellar age of 
45 Myr and a stellar mass of 3.7x 10 6 M Mpc -3 , about a factor of 3 lower than the minimum 
required to maintain an ionized IGM at z < 7. 



3.4. An Early Epoch of Reionization from Massive Stars ? 



The fundamental assumption here is that the galaxies seen at z ~ 6 are the remnants 
of star- formation responsible for reionization. It is plausible that the z ~ 6 galaxies that are 
seen in deep surveys are actually comprised of Population II stars with a Salpeter IMF which 
are responsible for a second late epoch of reionization. In that scenario, the earlier epochs 
would be entirely due to Population III stars which do not contribute to the ultraviolet and 
visi ble light luminosity density at z ~ 6 , a scenario that has been considered previously 



(e.g 



Cen 1120031 ; iFurlanetto fc Loeb 1 120051 ). These stars, which are more more massive than 
8M Q evolve into black holes within 30 Myr, i.e. by z ~ 6. The number of ionizing photons 
required to keep the IGM ionized between 15 < z < 6 is ~11 photons/baryon. An IMF 
extending between 10 — 100 M , with a Salpeter slope of 2.3 provides a total of 1.38 photons 
per baryon for a stellar mass density of 2.5xl0 6 M Q Mpc" 3 . The stellar mass density in 
massive Population III stars must therefore be 2xl0 7 M Mpc -3 . Thus, if reionization was 
initiated by massive stars which evolve into neutron stars and black holes by z ~ 6, there 
must be as much mass density in these remnants as in the stars that we detect in the galaxies. 

If in a contrived scenario, the initial z ~ 13 epoch of re ionization from Population 
III stars was relatively brief, lasting Az = 2 (e.g. ICen 1 120031 ). it would require that the 
stellar mass density in the initial burst was ~1O 6 M Mpc -3 . This corresponds to 10% 
of baryons that are in stars at z ~ 6. Since the lifetime of 10 M stars are ~30 Myr 
while the interval between 13 < z < 15 spans 60 Myr, the first epoch of reionization must 
be relatively inhomogeneous and will depend on the exact epoch at which star-formation 
in dark matter halos was initiated. Furthermore, the truncation of Population III star- 
formation th rough feedback process es argues against them being significant contributors to 
reionization (IGreif fc Bromm 1 120061 ) . Measuring the distribution of Stromgren sphere size 
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using high spatial resolution HI observations will reveal the true nature of the reionization 
history at z > 10. 

It should be noted that Population III star formation could potentially extend down to 
lower redshifts, depending on the effect of feedback on the metallicity of the star-forming 
environments. The rates of star-formation in such stars are however thought to be 3 x 10~ 4 
of the Population II star-formation rate and are inconse quential to the ultraviolet luminosity 
densi ty or the co-moving star-formation rate density (ITornatore et al. 1 120071 ; iBrook et al. 
2007h . 



4. Conclusions 

I have utilized existing data from deep Hubble and Spitzer surveys to provide a measure 
of the co- moving luminosity density at ultraviolet and optical wavelengths at z ~ 6. In 
particular, I provide a measure of the contribution to the optical luminosity density from 
faint galaxies which are below the Spitzer detection limit. Even after accounting for faint 
galaxies, the optical luminosity density is a factor of 2 — 3 below the ultraviolet luminosity 
density. I fit the resultant luminosity density estimates with a stellar population synthesis 
model to determine the maximal age and stellar mass density at z ~ 6. Assuming a Salpeter 
IMF, I find that the stellar mass density at z ~ 6 is 1.6 xlO 7 M Mpc~ 3 and the maximal age 
of a single stellar population is <100 Myr. By comparing the number of ionizing photons 
per baryon produced over the age of the starburst with the number of ionizing photons per 
baryon required to keep the IGM ionized, I find that reionization must have been a brief 
inhomogeneous process lasting <100 Myr and must have been completed as late as z < 7 if 
the stellar IMF had a Salpeter slope. 

Motivated by WMAP results which suggest an early epoch of reionization, I investigate 
the form of the stellar IMF if reionization was a single continuous process at redshifts higher 
than 6. If the past history of star-formation in z ~ 6 galaxies was responsible for reionization, 
I find that the the slope of the stellar IMF has to be non-Salpeter with a slope of a = —1.65 
if Zreion—9 an d at = —1.5 if z reion =ll for an IMF extending up to 200 M . However, the 
exact slope is sensitive to the ratio of the clumping factor to escape fraction, the metallicity 
of the stars and the value of the visible luminosity density at z ~ 6. 

On the other hand, the IGM could have been ionized between 6 < z < 15 by a population 
of massive stars between 10 — 100 M which evolve into stellar remnants like black holes by 
z ~ 6. Such a population would not be constrained by the rest-frame visible and ultraviolet 
luminosity density at z ~ 6. The resultant comoving mass density of remnants at z ~ 6 
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would be 2 x 10 7 M Mpc 3 and would therefore be comparable to the stellar mass density 
that we see in z ~ 6 galaxies. 

Final l y, I c onsid er the evolving stella r initial mass function that has been suggested 
by iDave I (120081 ) and Ivan Dokkum I (120081 ) to account for the discrepancy in specific star- 
formation rate estimates and colors of early-type galaxies at lower redshifts respectively 
This IMF if evolved out to z ~ 4, would result in a bottom-light IMF with a non-Salpeter 
slope at M < 12 M . Such an IMF could produce sufficient ionizing photons to account for 
late reionization of the IGM at z < 7 only if the visible light luminosity density was at the 
high end of our assumed range. 
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Table 1. Stellar Mass Density at z ~ 6 for Different IMFs 



Mass Range 


IMF Slope 


Stellar Mass Density 


Stellar Agc a 


Ionizing photons per baryon b 


M© 


a 


10 6 M Q Mpc~ 3 


Myr 




0.1 - 200 


-2.3 


16.1 


95 


0.91 


1 - 200 


-2.3 


7.3 


96 


2.06 


2 - 200 


-2.3 


6.1 


100 


2.69 


5 - 200 


-2.3 


5.1 


68 


3.97 


10 - 200 


-2.3 


2.5 


25 


5.50 


1 - 200 


-2.0 


15.1 


93 


3.36 


1 - 200 


-1.7 


9.5 


94 


4.58 


0.5 - 200 


-2.0 


10.7 


94 


2.96 


0.5 - 200 


-1.7 


15.6 


92 


4.36 


1-50 


-2.0 


7.0 


95 


1.88 


1-50 


-1.7 


8.6 


93 


2.77 


1-50 


-1.4 


11.6 


91 


3.72 


1-50 


-1.0 


19.0 


88 


4.88 



a This is only a maximal age for a single stellar population assuming all the stars formed in a single 
instantaneous burst and were then allowed to passively evolve. The exact age is sensitive to the e-folding 
timescale of a starburst and cannot be derived with the data at hand. 

b This is the total number of ionizing photons produced over the duration of the starburst per baryon 
per 10 6 Mq of stars formed. This needs to be multiplied by the escape fraction before a comparison 
with the recombination rate is made. 
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Table 2. Two Component Stellar Mass Density at z ~ 6 for Different IMFs 



Mass Range 
M 


IMF Slope 
a 


Maximal Stellar Mass Density 
10 6 M Q Mpc- 3 


0.1 - 200 


-2.3 


46 


1 - 200 


-2.3 


21 


2 - 200 


-2.3 


47 


5 - 200 


-2.3 




10 - 200 


-2.3 




1 - 200 


-2.0 


36 


1 - 200 


-1.7 


73 


0.5 - 200 


-2.0 


39 


0.5 - 200 


-1.7 


74 


1-50 


-2.0 


26 


1-50 


-1.7 


42 


1-50 


-1.4 


73 


1-50 


-1.0 


176 



c Maximal stellar mass density estimates at z ~ 6 derived using 
a two component stellar population. The young stellar component 
has an age of 10 Myr and dominates the ultraviolet luminosity 
density. The older stellar component has an age of 770 Myr and 
dominates the stellar mass density. See text for details. 
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Rest Wovelength (/^m) Rest Wovelength (//m) 

Fig. 1. — Left Panel: Solid circles show the rest-frame 1500A and 5500A luminosity density at 
z ~ 6 after appropriate completeness corrections for faint galaxies have been made. The faint 
galaxies have been assumed to have visible/ultraviolet luminosity ratios similar to that of 
z ~ 6 galaxies detected in Spitzer/ GOODS observations (the "high-V" case). The solid black 
line shows the age and stellar mass density of the template SED which fits the luminosity 
densities for a Salpeter stellar IMF with no dust extinction. The age indicates that if the 
stellar mass function had been Salpeter, the stellar population that is seen in z ~ 6 galaxies 
must have formed as late as 100 Myr earlier i.e. at z ~ 6.5. If faint galaxies have lower visible 
to ultraviolet flux ratios compared to bright galaxies, the V— band luminosity density would 
be decreased by 30% (the "low-V" case). The purple line is the spectral energy distribution 
which fits the luminosity density values in this latter scenario. Right Panel: Maximal fits to 
the "high-V" luminosity density at z ~ 6 assuming two stellar components and a Salpeter 
IMF. The UV luminosity density is fit by a 10 Myr old burst while a 770 Myr old, high 
mass to light ratio stellar component, is the dominant contributor to the visible luminosity 
density. 
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Fig. 2. — Left Panel: Comp arison between the clumping factors of neutral and ionized gas 
from the iTrac fc Cen I (120071 ) simulations for their reionization history. The blue dotted line 
shows the clumping of neutral gas averaged over the simulation volume. The green dashed 
line is the clumping of ionized gas averaged over the simulation volume. The solid black line 
is the clumping of ionized gas after star-forming halos have been excluded. The dot-dash red 
line is the clumping of ionized gas after all halos are excluded. Since, the dumpiness of gas 
in the star-forming halos is already factored into the escape fraction, we use the solid black 
line as the best estimate of the clumping factor. Right Panel: The plot shows the number 
of ionizing photons per baryon required to ionize the IGM and account for recombinations. 
The solid black line is the number of ionizing ph otons per baryon as suming the reionization 
history and redshift dependent clumping factor of ITrac &: Cen I (120071 ) after source halos have 
been excluded (solid black line in the left panel). The triple dot dashed line assumes that 
reionization was complete by z ~ 12 and the clumping factor and escape fraction remained 
constant at lower redshifts. The solid red line assumes the same reionization history but 
uses a constant clumping factor of 30. The dot-dash blue line assumes a constant clumping 
fa ctor of 10. The d ashed line is the number of photons derived using the parameteric solution 
of iMadau et al. I (119991 ) for c omplete hydrogen and helium reionization and the redshift 
dependent clumping factor of ITrac fc Cen I (120071 ). The reionization history and properties 
of the IGM are clearly important in determining the number of ionizing photons required. 
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Fig. 3. — Left Panel: The solid black line shows the minimum number of ionizing photons 
per baryon required to ionize the IGM and account for recombinations between the plotted 
redshift and z ~ 6 (Same as Figure 2). Overplotted as solid lines of different colors are the 
number of ionizing photons per baryon produced for a 0.02 Z Salpeter IMF (dN/dMocM -2 3 ) 
which reproduces the observed ultraviolet and optical luminosity density at z ~ 6. The 
different colors correspond to different low mass cutoffs; 0.1, 1, 2, 5, 10 M . The stellar mass 
density at z ~ 6 corresponding to these different IMFs are shown in Table 1. Reionization 
would be complete from the time at which the colored lines cross the solid black line. For 
the assumed clumping factor and escape fraction, a Salpeter IMF, even with a low mass 
cutoff, would maintain the ionized state of the IGM only at z < 7. Right Panel: The solid 
green lines show the number of ionizing photons for a 0.02 Z non-Salpeter IMF, with stellar 
masses spanning 1-200 M and different values of a, the IMF slope. The dashed blue line 
shows the number of ionizing photons for a 0.4 Z model template with the same range of 
slopes. Clearly, metallicity is an important parameter in constraining the IMF. 
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Fig. 4. — Dependence between stellar IMF and number of i onizing photons pro duced, as 
a function of IMF parameters for the rei onization history of iTrac fc Cen I (120071 ) . For the 
"late" reionization history considered by ITrac fc Cen I (120071 ). either the ratio of clumping 
factor to escape fraction must be much lower than the best estimates or the stellar IMF 
needs to have a slope index of a ~ —2.1 if the IMF extends up to 200 M . If the IMF has 
a high mass cutoff at ~50M Q , the slope needs to be a ~ —1.65. In either scenario, this is 
substantially flatter than a Salpeter slope (a—— 2.3). 
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Fig. 5. — Dependence between redshift of reionization and slope of the stellar initial mass 
function for the "high-V" case. The solid lines show the case where the IMF extends between 
1 — 200 M© while the dashed lines have the IMF extending between 1 — 50 M Q . Four different 
values of C/f esc are shown. The solid green lines are our best esti mate which assume s the 
redshift dependent clumping factor derived by the simulations of iTrac fc Cen I (120071 ) and 
an escape fraction of 0.1. The blue lines assume a redshift independent clumping factor of 
10 and an escape fraction of 0.15. The red lines assume a constant clumping factor of 30 
and escape fraction of 0.1 while the cyan lines are an extreme case assuming a clumping 
factor of 10 and an escape fraction of 0.5. For our best estimate of the clumping factor and 
escape fraction, the IMF must have a slope dN/dMoc M -L65 for reionization to be complete 
by z ~ 9 and dN/dMoc M -L5 for reionization to be complete by z ~ 11. 
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Fig. 6. — Similar to Figure 5 for the "low-V" case where the V-band luminosity density is 
30% lower as shown in Figure 1. Since a lower V-band luminosity density implies that fewer 
stars have been produced, the slope of the stellar IMF has to be even flatter resulting in 
even more high mass stars compared to Figure 5. 



